Are Proxima and Alpha Centauri Gravitationally Bound? 



O 
O 

(N 



O 



Jeremy G. Wertheimer 
Gregory Laughlin 

\J CO /Lick Observatory, University of California, Santa Cruz, CA 95064 

ABSTRACT 



Using the most recent kinematic and radial velocity data in the literature, we 
calculate the binding energy of Proxima Centauri relative to the center of mass 
of the a Centauri system. When we adopt the centroids of the observed data, 
we find that the three stars constitute a bound system, albeit with a semi-major 
axis that is on order the same size as a Centauri AB's Hill radius in the galactic 
potential. We carry out a Monte Carlo simulation under the assumption that 
the errors in the observed quantities are uncorrelated. In this simulation, 44% of 
the trial systems are bound, and systems on the 1-3 a tail of the radial velocity 
distribution can have Proxima currently located near the apastron position of 
its orbit. Our analysis shows that a further, very significant improvement in 



Oh' the characterization of the system can be gained by obtaining a more accurate 

O ' measurement of the radial velocity of Proxima Centauri. 

CO 

Subject headings: stars: individual (Proxima Centauri, a Centauri) 

1. Introduction 

In 1839, the Scottish astronomer Thomas Henderson learned of Bessel's successful paral- 
lax measurement of 0.314"±0.020" for 61 Cygni (Bessel 1838). Bessel's result gave Henderson 
the confidence to report his own measurement of a 1.16" ± 0.11" parallax for the a Centauri 
AB binary pair (Henderson 1839) 1 . Henderson's distance determination of 0.86 pc for a 
Centauri is reasonably close to the modern estimate of 1.35 pc. a Centauri remained the 
closest known stellar system until the measurement of Proxima Centauri's parallax (Voute 
1917). 



1 The Henderson (1839) reference is not returned by an ADS reference search on Au- 
thor=Henderson. However, a scanned copy of Henderson's 1839 MNRAS notice is available at: 
http://articles.adsabs.harvard.edU//full/seri/MNRAS/0004//0000168.000.html 
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Proxima and a Centauri exert a special fascination by virtue of the fact that they are 
the Sun's closest neighbors. Their hold on the imagination is further strengthened by the 
remarkable similarities of the masses and ages of a Cen A and B to the Sun. Proxima lies a 
mere 15, 000 ± 700AU from the a Centauri binary pair, and it has a small relative velocity, 
AV = 0.53±0.14km s _1 with respect to a Cen. The likelihood of such a stellar configuration 
occurring purely by chance is less than 10~ 6 (Matthews & Gilmore 1993), and based on this 
incredibly improbable arrangement, it has been suspected that the stars constitute a bound 
triple system ever since Proxima's discovery (Innes 1915). 

The dynamics of Proxima and a Centauri have been most recently examined by Matthews 
& Gilmore (1993) and by Anosova, Orlov, & Pavlova (1994). Both groups of authors showed 
that the (then-current) kinematic observations of the stars implied a positive total energy 
for the three-star system. Matthews & Gilmore (1993) sought to explain this somewhat 
disturbing result by hypothesizing that the measured masses of a Centauri AB were too 
low. They argued that a 3 — a increase in the dynamical mass of the system would result 
in a bound configuration. Anosova, Orlov, & Pavlova (1994) suggested that the three stars 
are members of a "stellar moving group", thus rendering Proxima's apparently hyperbolic 
trajectory less of an unusual occurrence. 

In addition to the reasons presented by Matthews & Gilmore (1993) for the importance 
of understanding the dynamical condition of the a Cen system, Proxima could also play a 
role in volatile enrichment of any terrestrial planets that might be orbiting either member 
of the central pair. A study by Wiegert & Holman (1997) showed that terrestrial planets 
are dynamically stable when placed within 4 AU from either star of the a Centauri binary 
pair. Furthermore, accretion calculations by Lissauer et al. (2004) suggest that terrestrial 
planet formation could have readily occurred within the a-Centauri system. One might 
therefore wonder about the habitability of putative terrestrial planets in the system. A 
possible concern with respect to habitability arises because any planets orbiting the a- 
Centauri binary may be depleted in volatiles. If Proxima were bound to the system during 
its formation stages, then it may have gravitationally stirred the circumbinary planetesimal 
disk of the a Centauri system, thereby increasing the delivery of volatile-rich material to the 
dry inner regions. 
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2. Data and Method 

The Hipparcos satellite returned extremely accurate kinematic information on all three 
stars of the a Cen system (ESA 1997). We can combine these data with mass and radial 
velocity information to formulate a complete dynamical picture of the system (see Table 1). 

The kinematic data in Table 1 are consistent with a broad spectrum of dynamical 
configurations for the three stars. We wish to determine what fraction of these allowed 
configurations correspond to bound systems. In order to do this, we perform an N=10000 
Monte Carlo simulation by assuming that each observed parameter is independent and varies 
with a normal distribution implied by its 1-a error. The only exception to the independence 
of the observations is that the parallaxes of a Cen A and B are set equal to each other within 
each realization, but still vary in a normal distribution from one realization to the next. We 
use observations of both a Cen A and B and the masses of each to calculate the properties 
of the center of mass of the A/B system. We calculate the gravitational binding energy E tot 
between the center of mass of the a Cen A/B binary and Proxima Centauri. We did this by 
converting to Cartesian coordinates and using cgs units with the simple formula 

E to t = Ta/b + T C + Ua/b/c (1) 

where T A /b and T c are the kinetic energies of the a Cen A/B center of mass and Proxima 
Centauri respectively and Ua/b/c is the gravitational potential energy of the two point mass 
system. 

3. Results 

We first calculated the properties of the a Cen A/B/C system using our method, but 
with the old observational data used by Anosova, Orlov, & Pavlova (1994). Following 
Anosova, Orlov, & Pavlova we also restricted our errors in the N=1000 Monte Carlo simu- 
lation to within their 2-a values. We find the A/B/C system is unbound with a probability 
of P=1.0 and that using the non-cgs (simulation) units from Anosova, Orlov, & Pavlova, 
the energy of the system was 5 ± 5 which was consistent with their results of P=1.0 and 
energy=6 ± 6. Using more accurate data from Hipparcos (ESA 1997), Pourbaix et al. 
(2002), and Queloz, D. (2004), not limiting the Monte Carlo errors to 2-a and using cgs 
units in an N=10000 simulation, we find that the binding energy of the A/B/C system is 
now 5 * 10 40 ± 17* 10 40 Joules with an unbound probability of P=0.55. With modern data it 
is now clear that almost half the realizations result in a bound triple system. Using the cen- 
troid of each observation we can compute an estimate of the observed orbital characteristics 
of Proxima Centauri with respect to the center of mass of the A/B binary (see Table 2). 
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For the case of a highly eccentric Keplerian orbit, Proxima Centauri would spend most 
of its time near apastron. The current distance of Proxima Centauri from a Cen A/B is 
15000 ±700 AU, thus the estimated semi-major axis of 272212 AU is almost certainly far too 
large. The velocity of Proxima Centauri relative to a Cen A/B in the Monte Carlo simulation 
is 0.53 ± 0.14 km s -1 . Systems with a lower relative velocity generally result in a bound 
orbit while higher velocities tend to result in a hyperbolic orbit (see Figure 1). According to 
observational data, Proxima Centauri is right on the cusp of being unbound, but we expect 
to find Proxima Centauri near apastron, which requires a slightly smaller relative velocity. 
Only P=0.025 of the runs had a low enough relative velocity to situate Proxima Centauri 
fairly near apastron. In Figure 2, we can see that the radial velocity of Proxima Centauri 
largely determines the energy of the resulting system and that a change of only 1 to 3-a in 
that radial velocity can result in a bound system currently near apastron. 



4. Conclusion 

The availability of Hipparcos data has provided us with the ability to implement a 
significant improvement over previous studies of the Alpha Centauri system. Our results 
indicate that it is quite likely that Proxima Centauri is gravitationally bound to the Alpha 
Centauri AB pair, thus suggesting that they formed together within the same birth aggregate, 
and that they three stars have the same ages and metallicities. As future observations bring 
increased accuracy to the kinematic measurements, it will likely become more obvious that 
Proxima Centauri is bound to the a Cen A/B binary and that Proxima Centauri is currently 
near the apastron of an eccentric orbit, (see Figure 3). Based on the expectation that the 
actual system is both bound and near apastron, we predict that improved measurements of 
Proxima Centauri's absolute radial velocity will yield a value of -22.3 km s _1 < v r < —22.0 
km s _1 . 

Acknowledgment: This work was supported by the US National Science Foundation 
CAREER Program under grant No. 0449986 to Greg Laughlin. 
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Fig. 1. — Histogram of the number of bound and unbound orbits of Proxima Centauri 
around a Cen A/B in the Monte Carlo simulation. The solid histogram represents the 
bound orbits while the dashed histogram represents the unbound orbits. The dotted vertical 
line represents the centroid relative velocity, which yields a barely bound orbit. 
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Fig. 2. — Distribution of total system energy versus radial velocity of Proxima Centauri 
in each Monte Carlo realization. The asterisks are orbits where Proxima Centauri is near 
apastron. The diamond at the center is the centroid value from the observations. The 
vertical dashed lines represent the 1, 2, and 3-sigma deviations of radial velocity. 
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Fig. 3. — Two simulated orbits of Proxima Centauri around a Centauri in the plane of 
the sky as seen from Earth. When Proxima Centauri has a radial velocity 1-a below the 
measured value, it results in a larger orbit while a 2-cr radial velocity results in a smaller 
orbit. Notice that in the 2-cr orbit, Proxima is closer to apastron which we would expect. 
The and 3-cr orbits are too large to be displayed on this scale. 
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Table 1. Observed Parameters and Masses for Proxima and a Cen A/B. 



Parameter 


Proxima 




a Cen A/B a 


ref 


a (RA) 


217°.44894751 ± 1.31 mas 


219° 


.91753275 ± 51 mas 


1 


S (Dec) 


-62°.68135207± 1.51 mas 


-60° 


.83712790 ± 35 mas 


1 


parallax (mas) 


772.33 ± 2.42 




742.12 ±1.4 


1 


^ Q (mas yr" 1 ) 


-3775.64 ± 1.52 




-3642.53 ± 12 


1 


^ (mas yr" 1 ) 


768.16 ± 1.82 




697.25 ± 9 


1 


V R (km s- 1 ) 


-21.8 ±0.2 




-22.445 ± 0.0024 


2,3 


mass (M B ) 


0.107 ±0.0214 




2.039 ±0.009 


4, 3 



a Errors in the a Cen A/B column are from the simulation (see text). 



References. - (1) ESA 1997; (2) Queloz, D., private communication 2004; 
(3) Pourbaix et al. 2002; (4) Henry 1999. 
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Table 2: Calculated Orbital Elements of Proxima Centauri around a Cen A/B. 

a (AU) Ellipticity Inclination (°) u (°) a (°) b Mean Anomaly (°) 
272212.148 0.985 150.9 86.87 197.81 359.53 



a argument of periastron 

b position angle of the ascending node 

Note. — These values are from the centroids of the observations and not our best guess from the Monte 
Carlo simulation. 



